Order by Disorder in Spin-Orbit Coupled Bose-Einstein Condensates 
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Motivated by recent experiments, we investigate the system of isotropically-interacting bosons 
with Rashba spin-orbit coupling. At the non-interacting level, there is a macroscopic ground-state 
degeneracy due to the many ways bosons can occupy the Rashba spectrum. Interactions treated at 
the mean-field level restrict the possible ground-state configurations, but there remains an accidental 
degeneracy not corresponding to any symmetry of the Hamihonian, indicating the importance of 
fiuctuations. By finding analytical expressions for the collective excitations in the long- wavelength 
limit and through numerical solution of the full Bogoliubov- de Gennes equations, we show that the 
system condenses into a single momentum state of the Rashba spectrum via the mechanism of order 
by disorder. We show that in 3D the quantum depletion for this system is small, while the thermal 
depletion has an infrared logarithmic divergence, which is removed for finite-size systems. In 2D, 
on the other hand, thermal fluctuations destabilize the system. 



I. INTRODUCTION AND OVERVIEW 

Multicomponent condensates of ultracold atoms offer 
rich physical systems due to the interplay between su- 
perfluidity and internal degrees of freedom fl] . Recently, 
through the use of synthetic gauge fields, two-component 
bosons with spin-orbit (SO) coupling have been engi- 
neered in the ultracold laboratory 2J. SO coupling in 
solid-state materials has a long history and is responsi- 
ble for a variety of interesting physical effects, with no- 
table examples including the spin Hall effect |3| and topo- 
logical insulators 3]- In addition, SO-coupled materials 
have diverse applications including spintronics [5]. The 
newer bosonic counterpart of SO-coupled systems using 
ultracold atoms have no analog in solid-state systems and 
are thus expected to exhibit genuinely new physics. SO- 
coupled cold atomic systems have also received consider- 
able recent theoretical attention -121). investigating 
topics such as spin-striped states [3 [T^ [T^, fragmen- 
tation [HI [19], and the realization of Majorana fermions 

[7|[I11[IS]. 

Recent experiments [2] have realized a special combi- 
nation of Rashba [22] and Dresselhaus [23] SO coupling 
in ultracold atoms. There are also promising proposals 
to realize more general non-Abelian gauge fields like pure 
Rashba (c.f. [24]), or even SU(A^) gauge fields that pro- 
vide a toolbox for topological insulators [25j. The con- 
ceptually simple system of a Rashba SO-coupled Bose- 
Einstein Condensate with isotropic interactions (RBEC) 
has surprisingly rich physics. The non-interacting system 
has a macroscopic ground state degeneracy as shown in 
Fig. [T] Interactions at the mean-field level partially re- 
move this degeneracy, but there remains an 'accidental' 
degeneracy not corresponding to any underlying symme- 
try of the system. Specifically, mean-field theory predicts 
a superposition of condensates of opposite momenta with 
their relative amplitudes and phases unspecified. 

In this work we show how fluctuations remove this ac- 




FIG. 1. (Color online) The non-interacting (Rashba) energy 
spectrum of the Hamiltonian Eq. ([T| with fc^ = 0. The red 
circle indicates the degenerate lowest-energy single-particle 
states. 



cidental degeneracy and select a unique ground state (up 
to overall symmetries) through the mechanism of 'order 
by disorder' [55]. Although the phenomenon of order by 
disorder has been theoretically accepted and discussed 
within the context of classical spin models '77 , quan- 
tum magnetism [351 ^'^'^ ultracold atoms [30 - 33 , ex- 
perimental demonstrations are, at best, scarce [33]. In 
contrast to the original proposal [H], the degeneracy 
lifting we find is primarily quantum driven. We deter- 
mine the fluctuation spectrum by numerically solving the 
coupled Bogoliubov-de Gennes equations. The resulting 
modes are integrated over to obtain the free energy as a 
function of the relative condensate weights and temper- 
ature. With this we show that fluctuations select a state 
with all bosons condensing into a single momentum state 
in the Rashba spectrum. We estimate the energy split- 
ting per particle due to fluctuations for typical experi- 
mental parameters to be on the order of 100 pK. While 
this splitting is smaller than typical condensate tempera- 
tures, it is the total energy which determines the ground 
state, so this effect should be readily observable provided 
the RBEC model can be realized. 
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II. DEFINITION OF HAMILTONIAN AND 
MEAN-FIELD GROUND STATES 

The Hamihonian describing non-interacting bosons in 
3D with SO couphng reads 

5£o = / dr*\r) ( f- - ■ p) *(r), (1) 

J \ Zm m J 

where '4'(r) = (4'^(r), ^'^(r))-^ is a two-component 
bosonic field operator, p is the momentum operator, 
Q is the magnitude of the SO couphng, and (t± is a 
vector composed of Pauh matrices as (t± = (cTj., ctj,, 0) 
(we set h — 1). The SO couphng in Eq. ([T]) is equiv- 
alent to the Rashba form j22j through a 90° spin ro- 
tation. The single-particle eigenstates of Eq. ^ have 
spins pointing either parallel or antiparallel to their mo- 
menta in the xy plane, and up to a constant have energies 

= 2k((fci±Q)' + fc')' ^here = {k^,ky,Q). 
Clearly, there is a ring in momentum space of degenerate 
lowest-energy states with = and |k_L| — Q (Fig. [T]). 
Correspondingly, there is a macroscopic number of ways 
TV non-interacting bosons can occupy this manifold of 
states. 

For the interacting portion of the Hamiltonian we take 
the simplest SU(2) invariant form 

JCi„t = J dr(|[p(r)]^-Mp(r)), (2) 

where p{r) = ^\r)'^(r), /i is the chemical potential, 
and g — where a is an effective scattering length. 
At the mean-field level one replaces the operators by c- 
numbers ^{r) '4'(r). The states that minimize the 
kinetic energy, Eq. ([l]), are in general given by 

|kj_|=Q.fc,=0 V / 

where are arbitrary coefficients and tan(c/Jk) = k^/ky. 
Minimizing the interaction energy restricts the mean- 
field states of Eq. ^ to have a constant density, p(r) — 
^^(r)^(r) = pq. Placing this constraint on states in 
Eq. ([3|, one finds that \If'(r) can have at most two nonzero 
coefficients A^^ occurring at opposite momenta. This can 
be shown by setting each non-zero wavevector component 
of *^(r)*(r) to zero. Without loss of generality, we take 
the momenta to point along the x-axis and thereby ob- 
tain the state 

*(r) - [ae^^^ ( | ) + ^e^^^ (4) 

where |ap + |5p = 1. We can take a and b to be real 
and parametrized as a = cos(|) and b — sin(|) since 
changing the phases of a and b amounts to position dis- 
placements and overall phase shifts of \J'(r) in Eq. Q. 
The selection of (a, b) as a result of spin-symmetry break- 
ing interactions (which is resolved at the mean-field level) 
was worked out in |12| . In contrast, in this work there 
remains a degeneracy at the mean-field level. 



III. CALCULATION OF COLLECTIVE 
EXCITATIONS 

The degeneracy with respect to 6 is accidental, i.e. it 
does not correspond to any symmetry of the Hamilto- 
nian 'K ~ 'Kq + J{int ■ We thus expect quantum fluctua- 
tions about the mean- field state Eq. Q to remove this 
degeneracy and to select a unique ground state through 
the order-by-disorder mechanism. To this end, we write 
^(r) = ^'(r)+'0(r) and perform a Bogoliubov expansion 
of 3-C to quadratic order in '0(r). Up to a constant the 
interaction Hamiltonian becomes J^int — 2 I '^^ i^Pi^)] 

where 6p{r) ~ ^^{r)'ip{r) + 'ijji\r)^{r). It proves use- 
ful to transform to the variables x(r) — {x^{r),xi{r))'^ 
where x{r) = e^'^^ ^ e~'^"''^^ e^'^y ^ rp{r) for which the in- 
teraction Hamiltonian takes the simple form ^int — 

f /dr(xt(r)+xl(r))'. 

The full Bogoliubov Hamiltonian, ^Keog, can be writ- 
ten compactly in matrix form if we introduce the four- 
component vector 

*(r) = (x^(r),x^(r))^. (5) 

Then up to a constant independent of 9 we find that 
^Bog = 5/ dr*\r)M(r,p)*(r), where 

2 

M(r,p)=l®l|- + ^(l-|-*T,)®(l-Fcr,) (6) 

Qp 

(1 (E) <Ty cos(2(5x) — (Tz® (Te sin(2(5x)) . 

In this expression, all 0-dependence is included in erg = 
cos(6')cr^ -I- sm{9)crz and (8) is the Kronecker product. 
This Hamiltonian can be diagonalized using a symplec- 
tic transformation [35l [36] , which amounts to solving the 
Bogoliubov-de Gennes equations, 

f7M(r, p)vk„(r) = £k„Vk„(r), (7) 

for positive eigenvalues £kn, where r/ — cr^ (g) 1, and 
Vkn(r) is a four-component function. Because of the 
translational symmetries of M(r, p), the eigenvalues are 
labelled with band index n and momentum k in the 
Brillouin zone (BZ) defined as —00 < ky,kz < 00 
and —Q < kx < Q- As usual, the eigenvectors are 
normalized as (vk„|j7|vk'„') = / drV|{^„(r)r7Vk'„' (r) = 
'5kk'^riri'Sgn(£kn)- In practicc, Eq. Q is simplest to solve 
in momentum space. Since the momentum space repre- 
sentation of M(r, p) is an infinite matrix, in numerical 
calculations it must be truncated at high momentum and 
the eigenvalues of interest must be checked to be inde- 
pendent of the cutoff. 

In Fig. [2] we show the two gapless (Goldstone) modes 
for several values of 9, found numerically from Eq. ([7|. 

In experiments of £q = ^ ~ gpo, so we set these 
quantities to be equal. The dispersion is plotted along 
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FIG. 2. (Color online) (a) The dispersions for the density and 
spin Goldstone modes for three values of 9 for kx — = 0. 
(b) The average (arithmetic mean) of the density and spin 

modes. In both plots we have fixed gpo — 



ky smce, as can be seen from Eq. ([7|, the spectrum £k„ 
has no ^-dependence when ky = 0. We refer to the dis- 
persions as 'density' and 'spin' modes since they reduce 
to the known expressions £kd — \J £k(£k + Sgpo) and 
£ks = f k ill the limiting case of Q = 0, where ek = 5;^ is 
the free particle dispersion. One sees that upon increas- 
ing Q from zero to 7r/2, the spin mode decreases in energy 
while the density mode increases. This gives, in a sense, 
a competing effect in terms of which (a, 6) configuration 
is selected from fluctuations. Noting this, in the right 
panel we plot the average of the spin and density modes 
for each value of Q. One sees that the average is always 
lowest in energy for = 0. This indicates that the zero- 
point fluctuations from the Goldstone modes will select 
= state though things become more subtle for T > 0. 
Such a state, as can be seen from Eq. (|4]), corresponds 
to all bosons condensing into a single momentum state 
of the RBEC system. The order-by-disorder mechanism 
will be considered more quantitatively below. 

Analytical expressions for the dispersions and eigen- 
vectors of Eq. ([7]) can be found perturbatively in the 
long- wavelength limit ek ^ ^Qi9Pa- In this limit one 

finds £kd = y 25/90 (e^^^ -I- Ae^^ sin^(6')) and £k>, = 

\j et°'''+AEl"'"^in^"(9) ^'"'^ density and spin modes respec- 
tively where A = gpo/(4eQ -I- Sgpo) and k^z = ^Jkl -f k^. 
These agree well with the numerical results shown in 
Fig. [2] for small k except for two special cases which re- 
quire more careful analysis. In particular, for = 0, the 
density mode disperses quadratically along ky while for 
< < 7r/2 the spin mode disperses as ky along ky. 
Otherwise the density and spin mode have respectively 
linear and quadratic dispersions about their minima. It 
is interesting to compare these to the noninteracting en- 



ergies shown in Fig. [Tj which have quadratic and quartic 
dispersions about their minima. 

IV. QUANTUM AND THERMAL ORDER BY 
DISORDER 

Let us now consider the free energy due to quantum 
and thermal fluctuations described by Jfeog- It is useful 
to separate out the contribution from zero-point fluctua- 
tions and write F(0) = Fq{e) -f Ft(0) where 

F,{0) = \ E (8) 

keBZ,n 

Ft{e)^kBT J2 ln(l-e-^^-('')), (9) 

keBZ, ri 

and /3 = l/ksT is inverse temperature. Reminiscent 
of the zero-point photon contribution to the Casimir- 
Polder force [37], the purely quantum contribution Fg(9) 
written as it is diverges. This divergence can be reg- 
ularized by subtracting the free energy for a particular 
mean-fleld configuration which we take to have = 0: 
AFq{9) — Fq{6) — Fq{Q). This regularized expression 
converges 38 , and no renormalization of the effective 
range of interactions is needed. The zero-point contribu- 
tion to the free energy numerically computed as a func- 
tion of is shown in Fig. [sj^a) where the summation is 
performed over 26 bands (we emphasize that in order to 
obtain quantitatively correct results, including only the 
gapless modes is insufficient). One sees, indeed, that the 
= state has the lowest energy and at T = such a 
state is unambiguously selected. 

We now turn to the finite-temperature contribution to 
the free energy. Interestingly, one finds that the sign of 
the thermal contribution AFt{6) = Ft{6) — Ft{0) is neg- 
ative and opposite to that of AFq{6). Furthermore, the 
magnitude of the thermal contribution is always smaller 
than the contribution from zero-point fiuctuations, in 
contrast to typical situations where thermal fiuctuations 
enhance the degeneracy lifting and are larger in magni- 
tude for modest temperatures (see, e.g. [3Dj). Another 
instance of where quantum and thermal fluctuations se- 
lect different states is in Ref. [33]. The sign of AFt at 
low T can be understood by noting that the spin mode 
has the lowest energy for = 7r/2 (Fig. [2]). 

As seen in Fig. [sjjb), the magnitude of AFt{d) ap- 
proaches AFg{9) at high T, so that AF{e) = AF,(0) + 
AFt{0) = (T-i) ^ in this limit. This behavior can 
be understood through a high T expansion of the free 
energy 

Ft{0)^kBT E ln(/3£k„(0))-^ E £k„W. (10) 

keBZ,n keBZ,n 

As the second term cancels the quantum contribution. 
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(b) 




0,25 0.5 50 100 

FIG. 3. (Color online) (a) The zero-point contribution to the 
free energy AFq as a function of 9 for three different values 
of g. (b) The absolute value of the (negative) thermal free 
energy splitting between the = and = 71/2 configurations 
|A_Ft(7r/2)| as a function of temperature (solid line). This is 
seen to approach the quantum zero-point splitting AFq{n/2) 
at high temperatures (dashed line). In both panels the solid 
lines are fits to the numerically computed data points. 

we focus on the larger first term which can be written as 
ksT J2 ln(/?£k„(0)) = ^fcBrin|det(/3r7M(r,p)| 

keBZ.n 

keBZ'.n 

where Akn are the eigenvalues of M(r, p) and we have 
used |det(?7)| = 1. The second summation above is over 
the reduced Brillouin zone BZ' which is restricted to pos- 
itive values of k^. The eigenvalues Akn are independent 
of the condensate configuration given by 9. This can be 
seen by noting that the 0-dependence of M(r,p) can be 
removed through the unitary transformation M(r, p) — )■ 
U^M(r, p)U where U = l{l + (T^)^l + ^{l~(T^)(x)e'''''y 
[31]. Thus to this order we find that Ai^f (61) = -AFq{9). 
The next-order term in the high-temperature expansion 
has 1/T dependence which is evident in the numerical 
results shown in Fig. [sjjb) . 

V. CONDENSATE DEPLETION 

Having established using the Bogoliubov expansion 
that fluctuations select 6* = 0, we now investigate the 
self-consistency of this approach. This is determined by 
the depletion or the number of particles excited out of the 
condensate N^x considered as a fraction of the total parti- 
cle number N . Consistency of course requires that this be 
finite, while neglecting of terms beyond quadratic order 
in JfBog is quantitatively reliable only if A^ex ^ N. The 
quantum and thermal contributions to A'ex = Nq + Nt 
are, respectively, Ng = ^ I]keBZ,«('"kn|(l - ■^)|wkn> and 
= EkeBZ,«(^kn|wkn)/(£kn), where f{x) = (e*^^ - 
is the Bose-Einstein distribution function. The only 



possible divergences of these expressions are in the in- 
frared, and so can be studied analytically using the small- 
k expansion. 

In 3D at T = 0, Nq is finite and so can be sufhciently 
small provided weak-enough interactions. Our numerical 
results in fact demonstrate that the depletion is small 
even for moderately strong interactions, including in the 
region of experimental relevance. For T > 0, the ther- 
mal contribution is instead found to have a logarithmic 
divergence in 3D. This divergence (similar to that occur- 
ring in quasi-2D scalar condensates) is naturally removed 
for finite-sized systems and the condensate will thereby 
satisfy the stability criterion. For 2D condensates with 
isotropic SO coupling, the situation is different. Here, 
the quantum depletion again converges, but the thermal 
depletion diverges as 1/k for small k. Thus, at T > 
our theory is unstable in 2D, which is consistent with 
work suggesting fragmentation [8j . Our conclusions 
on the stability of the condensate are, remarkably, iden- 
tical to those based on the simple application of the Ein- 
stein criterion to the noninteracting spectrum. This is 
particularly surprising since the low-fc region is strongly 
modified by the interactions, giving quasiparticle modes 
that disperse with different powers than in the noninter- 
acting case. 



VI. EXPERIMENTAL FEASIBILITY 

We now comment on the experimental feasibility of 
observing order by disorder in RBEC. We first consider 
the magnitude of the degeneracy lifting. As a prototyp- 
ical example we take spin one ®^Rb. For a typical den- 
sity of po = 2 X lO^^cm^'^ and gpQ = we find (for 
the appropriate scattering lengths) that at zero temper- 
ature the free energy splitting per particle due to fluc- 
tuations is AF{TT/2)/kBN = 110 pK. One should note 
that this number should not be directly compared with 
the condensate temperature since the total energy de- 
termines the ground state. It is this energy which will 
determine experimental timescales for the relaxation to 
the ground state. Spin-one atoms also possess a spin- 
dependent interaction term which we have neglected. 
For ^^Rb, however, this spin-dependent interaction is 
relatively small (.5% that of the spin-independent) and 
as a result the degeneracy lifting from fluctuations is 
typically larger. Alternatively, one could use fermionic 
homonuclear molecules that have a singlet ground state. 
More importantly, schemes to create SO coupling in 
bosonic systems typically rely on utilizing dressed states 
[31ITH1I21] which can induce anisotropic interactions. The 
magnitude of such terms and their effect on the order-by- 
disorder mechanism will be left to future work when it 
becomes clear which of the several proposed schemes is 
most promising to realize Rashba coupling. 

Another entity of experimental relevance is the har- 
monic conflning potential. The results of the current 
manuscript will hold if the conditions for the local den- 
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sity approximation (LDA) are satisfied [30] • For our sys- 
tem this requires that the energy sphtting of the single- 
particle states as recently found in in pOl be small 
compared to the interaction energy. This energy split- 
ting becomes small for weak trapping and/or strong SO 
coupling, resulting in a large quasi-degenerate manifold 
of single-particle states. Conversely, in the weakly in- 
teracting limit (where LDA is inapplicable) recent work 
[inim] has shown the ground states of the RBEC system 
in a trap can form vortex lattices. 

VII. CONCLUSION 

In conclusion we have investigated the system of 
Rashba SO-coupled bosons with isotropic interactions 
(RBEC). In general bosons with SO couphng offer a gen- 
uinely new class of systems which has not been addressed 
in the vast solid state literature on spintronics [5^ . In par- 
ticular, we have established that fluctuations select the 
RBEC system to condense into a single momentum state. 
We have argued that such a configuration is stable in 3D 



but destabilized when T > in 2D. We expect bosons 
with Rashba SO coupling to be realized in the near future 
for which the predicted configuration will be observable 
in Stern-Gerlach experiments. In future studies it will 
be interesting to investigate more general combinations 
of Rashba and Dresselhaus SO coupling. Such systems 
also possess accidental mean-field degeneracies and thus 
fluctuations are expected to play an important role in 
determining their ground states. In addition it will be 
interesting to investigate RBEC systems in two dimen- 
sions. 
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